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1. SINGLE PHONON MODES 



1.1 Introduction 

Mechanical vibrations are described as a collection of normal modes. When these 
are quantized, we obtain phonons - very similar to the single harmonic oscillator. 
Phonons are thus quantum excitations. A single phonon can involve billions of 
atoms, i.e. all those forming the original normal mode. 

In bulk materials, one usually finds a whole spectrum of normal modes, char- 
acterized by a dispersion relation 



where k is a label for the normal mode (usually a <i-dimensional wave vector, were 
d is the dimension of the vibrating system), and uj is the angular frequency of the 
normal mode (two pi times the frequency). We will deal with bulk phonons in the 
second part of this Lecture. 

Sometimes certain normal modes of a vibrating system are particularly inter- 
esting, for example the fundamental (lowest) vibration mode in a small beam or a 
string, or a vibration mode in a molecule. Often such vibrations are observed by 
interaction with other quantum objects, for example phonons (light) or electrons. 

The interaction between electrons and phonons is, needless to say, of key im- 
portance for our understanding of a broad variety of physical phenomena such 
as superconductivity, thermal conductance, electrical resistance, and decoherence, 
just to name a few. The chapter 'electron-phonon interactions' thus forms part of 
any good textbook on, e.g., solid state physics (such as 'Ashcroft/Mermin' just for 
the sake of naming one...). 

The interactions between electrons and phonons are however also a paradigm 
for what one would call interactions between massive particles and gauge bosons 
in particle theory. I am interested in the 'cleanest' possible kind of such interac- 
tion, and that is why the subject of this first Lecture is the interaction between 
single electrons and single phonon modes. At first sight, this might perhaps look 
somewhat artificial but we will soon see that some of the most exciting, recent 
experiments in nanomechanics, molecular electronics, and mesoscopic physics are 
just about this idea. 



u = w(k) 
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1.2 A Simple Model Hamiltonian 

Rather than giving an overview of recent experiments with single phonon modes, 
I want to move directly to the definition of the simple model mentioned above. 
It describes electrons tunneling between a left and a right lead (electron reservoir 
7i ICS ) through a single molecular level e that is coupled to a single phonon mode 
with creation operator at. The total Hamiltonian (system plus coupling to leads) 
is 

H = H s + H v +H rcs (1.2.1) 

The system Hamiltonian Tis describes a single molecular level e (we will sometimes 
call it 'dot') linearly coupled to a single harmonic oscillator mode, 

H s = uj a J a + \u (a + a))n d + en d + ^-n d (n d - 1). (1-2.2) 

where n d is the electronic number operator. The term with U mimicks the Coulomb 
interaction between two electrons with spin up and sin down on the level. In 
the following, we assume that only one electron can occupy the molecular level, 
corresponding to the limit U — > oo. We also assume that the spin of the electron 
plays no role here - this is a little bit artificial as the electron spin does play 
a role usually. We could imagine, for example, that the electrons are all spin- 
polarized in a very strong magnetic field. The point is that we really want to 
focus on the electron-phonon interaction here and therefore do not wish to include 
any additional interaction effects in our mdel. In the following we will stick with 
spin-polarized electrons (e.g., all electrons with spin f) and omit the spin label. 
As we have only n d = or n d = 1 electrons on the dot, we shall drop the term 
Y n d( n <i — 1) from now on. The electron number operator on the dot is 

n d = Sd= |1)(1|, d = s = |0)(l|, (1.2.3) 

where d is the annihilation operator for the electron on the dot and we used an 
equivalent representation in Dirac bra-ket notation with |0) denoting the empty 
and |1) denoting the occupied dot state, writing d as a projector (for which we 
reserve the symbol s). 

The electron reservoir term is simple, 

k,a=L/R 

describing a gas of free spin-polarized electrons in the left and right reservoir. 
Finally, the tunnel coupling between the reservoir and the dot is 

H v = J2( V ^ c lj + H-c), 8 = \0)(1\. (1.2.5) 

ka 

Here and in the following, we always set h — 1. 
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1.2.1 System Hamiltonian 

We first focus on the system Hamiltonian. 

Harmonic Oscillator 

The vibrational part of the eigenstates for rt^ = are the usual harmonic oscillator 

states 

Ho\n) = nu \n) , H = uj a^a, (1.2.6) 

which are n-phonon states corresponding to a vibrational energy nu (we shifted 
the energy scale such that the \?vjJq from the ground state energy no longer ap- 
pears). 

Shifted Harmonic Oscillator 

For nd = 1, the vibrational part of Tis describes a shifted harmonic oscillator, 

H\ = ujqo) a + \uo G {a + a)) = ujQ^b — ujq\ 2 (1.2.7) 
b = a + A, b f = a^ + \, (1.2.8) 

where we introduced shifted ladder operators b and W (note that [b, b'] = 1 = 
[a, a']). The operators a and b are related by a unitary transformation, 

b = a + A = X x aX{, X x = e x(a ~ at \ (1.2.9) 

which is shown by a nested commutator expansion (see below). The operator 
X\ is called unitary displacement operator. We use it to unitarily transform the 
unshifted Hamiltonian Hq into the shifted Hamiltonian H\, 

X X H X{ = H X + oo X 2 = Uotfb. (1.2.10) 

The eigenstates of H\ are given by shifting those of Hq, 

H \n) = nu \n) ~> H x X x \n) = [nuj - Ld X 2 ]X x \n) , (1.2.11) 

which means that the eigenstates of H x are shifted oscillator states 

\n) x = X x \n) (1.2.12) 

with energies uujq — uqX 2 . In particular, the groundstate of H x is a coherent state 
in the basis of the unshifted oscillator with energy 

s = -uj \ 2 . (1.2.13) 
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EXERCISE: Why is the ground state energy of the shifted oscillator smaller than 
the ground state energy of the unshifted oscillator? Hint: draw a picture of the 
oscillator potential in the position representation. 

Coherent states |a)o are eigenstates of the phonon annihilation operator a with 
complex eigenvalues a, 

a\a)o = a\a)o, a. G C (1.2.14) 

and are usually labeled by Greek symbols. We have 

6|0) A = a\0) x = — A|0> A ~» |0) A = X A |0> = | - A) . (1.2.15) 

The ground state of H\ thus is a coherent state |a)o in the basis of the unshifted 
oscillator with a = — A. 

Diagonalized System Hamiltonian Tis 

The diagonalized system Hamiltonian Tis with two electronic states, empty |0) e 
and occupied |l) e , therefore has eigenstates 

\0n) = |0) e ® \n) , | In) = |l) e ® |n) A . (1.2.16) 
1.3 The Polaron Transformation 

The polaron transformation is a general method to treat problems with linear 
coupling between electrons and phonons, such as in our model 7i, Eq. (1.2.1). We 
discuss the method first for our system Hamiltonian Tis, Eq. (1.2.2), which, of 
course, will reproduce our previous results for the shifted oscillator. 

1.3.1 System Hamiltonian 

We introduce a unitary transformation for all operators O, including the Hamil- 
tonian, 

d^O = e s Oe~ s , S = n d X{a^-a). (1.3.1) 

The operator U = e s is unitary because U^ 1 = U\ i.e. UW = U^U = 1. Note 
that rid in the definition of S is the electron number operator and not just a number 
(0 or 1) which is important when it comes to transforming d and d* for the dot 
electrons. 

For explicit calculations, one has to use the nested commutator expansion 1 , 
6 = O + [S,0] + I [5, [S,0]] + [5, [5,0]]] + ... (1.3.2) 

1 The nested commutator expansion is proved by solving a first order differential equation for 
fix) = e Oe , using a power series in x and comparing coefficients. 
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We first transform the electron number operator n d , 

n d = n d + [S, n d ] + ... = n d , (1.3.3) 

because S and n d commute. This is not the case for the phonon annihilation 
operator a, 

a = a + [S, a] + -^[S, [S, a]] + .... = a + n d \[a) - a, a] + .... 

= a-n d X, (1.3.4) 

because [a*, a] = — 1 is a scalar and all the higher commutators disappear. Corre- 
spondingly, 

t _ „t _l rc „ti i J_rc re „t 



a' = a 



2! 

= a' — n d A, (1.3.5) 

as we expected by taking the canonical conjugate of a. 

Transforming now all operators individually, the polaron transformed Tis is 

H s = e s H s e~ s = cu e s a t e~ s e s ae~ s + \u e s (a + a t )e~ s e s n d e~ s + ee s n d e~ s 

= ujQa)a + \ujq(cl + c^)n d + en d 

= io ((a f - n d \) (a - rid A)) + Xu {a - n d \ + a f - n d X)n d + en d 

= ujqo) a + en d — Lu \ 2 n d . (1.3.6) 

The square of the number operator is just the number operator itself, n 2 d = n d . 
there is either one or no electron on the level (n d thus is a projection operator that 
projects onto the number of electrons on the level). We can thus write 

Hs = lu^ a + (e — LU \ 2 )n d . (1.3.7) 

The effect of the polaron transformation is thus twofold: 

• the transformed Hamiltonian Tis is diagonal in the (old, unshifted) phonon 
operators a and a*: the linear term with the coupling to n d has been trans- 
formed away. 

• the energy e has been shifted into a lower energy e — ujq\ 2 (the phonon 
frequency is always positive, u;o > 0). 

How do we now use the polaron transformation in oder to obtain the eigenvalues 
e n and eigenstates \n)\ of 7is? As Hs is diagonal in the old phonon operators, one 
obviously has 

H s \n) = (nu + e - u \ 2 )\n) , (1.3.8) 
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where |n)o are the usual eigenstates of the phonon number operator a^a. This 
is, however, still in the polaron-transformed frame. We obtain the eigenstates 
\n)x of the original Tig by back-transforming: write the eigenvalue equation for 
vibrational states of the occupied level rid = 1 as 

T-(-s\n)\ = e n \n)x,^Hse s \n)x = e n e s \n) x 
ujQ(vae \n)x = (e n — e + uo\ 2 )e \n)\ (1.3.9) 

Comparing with the usual, unshifted harmonic oscillator eigenvalue equation yields 

e S \ n )\ = \n)o, nu = e n — e + lu X 2 . (1.3.10) 

The eigenvalues e n and eigenstates \n)x of Tis are therefore 

\n) x = e- s \n) = e x{a - ar) \n) , e n = mj + e - u \ 2 . (1.3.11) 

1.3.2 Total Hamiltonian 

Electron transport through the system is introduced by connecting to left and 
right electron reservoirs in thermal equilibrium at chemical potential \xl (i^r) with 
positive source-drain voltage Vsd = — H>r > 0, inducing tunneling of electrons 
from the left to the right. We recall that the total Hamiltonian is Ti, = Hs + Ti-v + 
Tires, Eq. (1.2.1). We recall that the total Hamiltonian (system plus coupling to 
leads) is 

7~t = Hs + Hv + Tires, Tires == ^ £fc,aC fea ,Cfe Q (1.3.12) 

k,a=L/R 

Tiv = J2( V *<x c la § + H.c), s = |0)(l|. (1.3.13) 

ka 

We have to transform s = |0)(1|, and using the commutator 

[n d ,s] = |1)(1|0)(1| - |0)(1|1)(1| = -s, (1.3.14) 
we obtain from the nested commutators with S = n^A(a^ — a), 

s = e s se- s = s + [S,s] + ^[S,[S,s}} + .... 

1 2 

= s — A(a^ — a)s + — (X(a' — a)) s — ... = Xs 

X = e A(a " at) (1.3.15) 
The transformed Hamiltonian thus is 

H = u a^a + {e-u Q \ 2 )n d + Y^{V k J k JX + H.c.)+Tires- (1.3.16) 

ka 

We note that the electon reservoir Hamiltonian, 7i rGS , is left invariant under the 
polaron transformation. 
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1.3.3 Polaron-transformed Master Equation 

The Master equation is now derived in second order perturbation theory in the cou- 
pling to the leads. Performing the commutators and using the free time evolution 
of the electron reservoir operators, one finds 

jp(t) = - y, f dt'gut-t r ){msKmt')-~s(typit')m} 

ka=L/R-*° 

~ T fdt'hait'-t) {s\t)s(t')p(t') - S(t>)p(t')^(t)} 
ka J ° 

- E f t dt'g ka (t'-t){p(trs(trsKt)-~s^mmt')} 

, Jo 

ka 

- T fdt'g ka {t - {~p(t')~sKt')~s{t) - s(t)p(t')S\t')} 

ka J ° 

g ka (r) = \V ka \ 2 f a (e k )e^ T , g ka (r) = \V ka \ 2 [l - f a (e k )}e^\ (1.3.17) 

with the Fermi distributions f a (e ki ) = Tr rcs (R cl a c ka ) . Here, the time evolution 
of the polaron-transformed s-operators is 

s(t) = e inst se- inst . (1.3.18) 

The sums over k can be written as integrals, introducing the tunneling density of 
states v a (e) in lead a, 

/oo 
deu a (e)r (8)6^1 u a (e) ee £ \Vka\ 2 5(e - e k ). 
■°° k 

(1.3.19) 

The further calculation is now performed as follows: 

1. Take matrix elements of Eq. (1.3.17) by 'sandwiching' as 

(0,n|...|0,n'}, (l,7i|...|l,n'}, (1.3.20) 

where |0n) = |0) e ® \n)o etc. as in Eq. (1.2.16). In particular, the phonon 
numbers n refer to the un-shifted harmonic oscillator. 

EXERCISE: Check that there are no terms (0, n|...|l, n'), (1, n|...|0, n'). This 
means that there are no coherences between occupied and un-occupied states. 

2. Explicitly insert the time-dependence of the s(t') etc. 

3. Perform the Markov approximation (if required) in the interaction picture. 
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4. Transform back into the Schrodinger picture via 



(0,n\p(t)\0,n') 



(0,n\p(t)\0,n') 



1.3.21) 



d 

7ft 



(0,n\p(t)\0,n') = z(eon'-eon)<0,n|/>(t)|0,n') 

d 



+ e i{£ ^'- £0nt) — (0,n\p(t)\0,n'). (1.3.22) 

(JjL 

5. Check that in all exponents the time-dependence is always via t — t'. 
This derivation yields 

d ., 



dt 



— H £ 0n' — £Qt, 



(0,n|p(t)|0,n') 



E 



dev n (e) 



dt' x 



f a (e)(0,m'\p(t)\0,n')(n\X\m) (m\X^\m')e i( - £ - £lm+£ °™' ){t ~ t ' ) 
+ f a {e) (0, n\p{t)\0, m) (m\X\m) (m > \X^\n')e- i{£ ~ £ ^' +£0m){t ~ t ' ) 



f a (e)(l,m\p(t)\l,m')(n\X\m)(m'\X^\n')e 
f a (e)(l,m\p(t)\l,m'){n\X\m)(m'\X^\n') 



-i(e-e lm +e 0n )(t-t') 



'\J( £ - £ im'+ e On')( t - t ') 



1.3.23) 



and a corresponding equation for the (1, n\p(t)\l, n'). The Pranck-Condon Fac- 
tors 



(n\X\m) 



1.3.24) 



etc. are determined by the overlap of the vibrational wave functions before and 
after tunneling of one electron. They strongly determine the electronic transportin 
our model, in particular if the electron-phonon coupling is strong which means 
A ~ 1 or larger. 



1.4 Franck-Condon Factors 



1.4.1 Born-Oppenheimer Approximation in Molecular Physics 

(This is from my lectures notes 'Quantum Mechanics of Atoms and Molecules'). 

Molecules are system consisting of electrons and nuclei. This definition cov- 
ers the full range from rather simple molecules like H2 up to extremely complex 
situations with billions of nuclei, or in principle even solids or fluids although one 
usually thinks of something like a microscopic object. The question, of course, 
is what microscopic really means. In principle, one could have molecules with 
macroscopic large numbers (like 10 23 ) of electrons and nuclei. Would these behave 
as quantum or as classical objects? 
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Model Hamiltonian 

We start from a Hamiltonian describing a system composed of two sub-systems, 
electrons (e) and nuclei (n) 



where H en is the interaction between the two systems. Note that the splitting of 
the Hamiltonian H, is not unique: for example, Ti, n could just be the kinetic energy 
of the nuclei with their mutual interaction potential included into 7i en (as in the 
BO approximation). 

The set-up 7i = H e +TC n +li. cri is quite general and typical for so-called 'system- 
bath' theories where one would say the electrons are the 'system' and the nuclei are 
the 'bath' (or vice versa!). In the theory of molecules, however, things are a little 
bit more complicated as there is a back-action of from the electrons on the nuclei. 
This back-action is due to the electronic charge density acting as a potential for 
the nuclei. 

There is no a priori reason why the nuclei and the electronic system should 
not be treated on equal footing. However, the theory has a small parameter 



given by the ratio of electron mass m and a typical nuclear mass M 3> m, and the 
exponent 1/4 is introduced for convenience in the perturbation theory used by Born 
and Oppenheimer in their original paper. The smallness of this parameter makes 
it possible to use an approximation which is called the Born-Oppenheimer ap- 
proximation. 

We assume there is a position representation, where q = {x 1; x^} represents 
the positons of all electrons, X = {Xi,...,Xn} the positions of all nuclei, and 
correspondingly for the momenta p and P, 



Spin is not considered here. Also note that the interaction only depends on (q, X) 
and not on the momenta. 

The Born-Oppenheimer Approximation 

This is the central aproximation used in many calculations. We now try to attack 
the Schrodinger equation = £^ for the total system (electron plus nuclei). 

A first guess to solve the stationary Schrodinger equation Ti^ = £^ for the 
total system would be a separation ansatz 






H = H(q, p; X, P) = H e (q, p) + H n (X, P) + H en (q, X). 




H(q,p;X,P)*(q,X) 



ifj e {q)4> n {X) unsuccessful, 
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which does not work because the interaction TC en (q, X) depends on both q and X. 
As Ti, ea (q,X) depends on the positions of the nuclei X, let us try an ansatz 

V(q,X) = ip e (q, X)(j) n (X) successful (1.4.5) 

where now the electronic part depends on the nuclear coordinates X as well. This 
looks unsymmetric: why shouldn't one have ^(q,X) = ip e (q, X)(p n (q, X)l First, 
there is an asymmetry in the problem in the form of M > m, and ty(q,X) = 
ip e (q, X)(j) n (q, X) is no more better than ty(q,X) in the first place. 

The idea with writing ty(q,X) = ip e (q, X)<fi n (X) is that the electronic part 
ip e (q,X) already solves part of the problem, i.e. 

iH c (q, P ) + H cn (q,x)]M<i,x) = e(x)Mq,x), (1.4.6) 

an equation in which X, of course, appears as an external classical parameter that 
commutes with all other variables. Consequently, the eigenvalue E(X) has to 
depend on X. We thus obtain 

Hlpe4>n = [H e + H Q + H cn ] tp e (f> n 

= [H n + E{X)}^ n (?)=£VW>„ (1.4.7) 

where the last questionmark indicated what we would like to have! Since H, n and 
E(X) depend on the nuclear coordinates only, one would like to use an equation 
like ' 

[Hn + E{X)\ MX) = £0nPO, (1.4.8) 

because then we would have achieved our goal. However, the operator 7i n contains 
the nuclear momenta P which operate on the X in ifj e (q,X), i.e. 

Hlpe4>n = Ipe [H n + E(X)} (j) n + [H n j)e4>n ~ ^ e H n n ] 

= £'lp e (pn + [H n ljj e (j) n - lp e H n (p n ] . (1.4.9) 

This shows that we are almost there if it wasn't for the underlined term. One now 
tries to find arguments why this term can be neglected. If it can be neglected, 
then we have achieved the full solution of the Schrodinger equation by the two 
separate equations 

[H e {q,p) + H eQ (q,X)]ip e (q,X) = E(X)ip e (q, X) electronic part 

[H u + E{X)\ <f> n (X) = £<f> n (X) nuclear part. (1.4.10) 

These two equations Eq. (1.4.10) are the central equations of the Born-Oppenheimer 
approximation. Even without solving them, some quite interesting observations 
can already be made: 
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Figure 1.1: Franck-Condon-Principle. Left: Classical picture, right: quantum- 
mechanical picture. From Prof. Ed Castner's lecture http://rutchem.rutgers.edu/. 

• The electronic part is calculated as if the nuclei were at fixed positions X 
('clamped nuclei'). 

• The eigenvalue of the energy of the electronic part serves as a potential energy 
for the nuclei in the nuclear part of the equations. 

The Franck-Condon Principle 

Here, a good description is in Atkins/Friedman ch. 11.4. 

We consider a transition a — > a' between two molecular states due to, e.g., 
an optical excitation. For simplicity, we leave out the molecular rotations here 
and just discuss electronic and vibrational transitions. In a classical picture (with 
respect to the large mass nuclear motion), one considers the two effective potentials 
U a (r) and U a >(r) acting on the nuclei as given functions. One then argues that 
the electronic transition occurs so fast that the nuclear system has no time to 
react: before and after the transition, the nuclear coordinate X is the same. This, 
however, means that the distance \x'\ = \X — r a > \ from the equilibrium position r a > 
after the transition and the distance \x\ = \X — r a \ from the equilibrium position 
r a before the transition are not the same: when the nuclei are in equilibrium 
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before the transition (X = r a ,x = 0), their new coordinate x' relative to the new 
equilibrium r a > is x' = X — r a > = r a — r a > ^ after the transition. 

The optical transition depends on the total dipole moment operator which is 
a sum of electronic and nuclear dipole moment, 

d = -e^ qi + e^Z s X s = d e + d n . (1.4.11) 

i s 

The transition matrix element in Born-Oppenheimer approximation is (a ^ a') 
(aV|d e + d n M = j dqdXrAqX)K>y(X)(de + d Q )<P atV m a (qX) 



dqip* a ,(qX)d e ip a (qX) 



dXf a , :V ,(X)^ v (X)d n / dqr a ,(qX)MqX) 



dXfcyiX) J dqrAqX)d e ^ a {qX) 

~ (a'\d e \a)S(v, v'), S(v, v ') = (v'\v). 
Here it was assumed that the integral 

dqtp* a ,(qX)d e ip a (qX) w (a'|d e |a) 



(f>a,v( X ) + 

(1.4.12) 
(1.4.13) 



does not depend on the nuclear coordinates X. 

The transition between two electronic levels a and a' is therefore determined 
by the dipole matrix element (a'|d e |a) and the Pranck-Condon factors S(v,v'), 
which are the overlap integrals of the corresponding vibronic states. As these 
states belong to different electronic states a and a', the overlaps are not zero, and 
there is also no selection rule for Av. 

1.4.2 Example: Linear Oscillator Coupling 

In this model, the Franck-Condon factors can be calculated analytically: they are 
given by the phonon number matrix elements of the displacement operator 



D(a) = e a{a] ~ a) 



1.4.14) 



and can be obtained using |n) = (l/^/n\)(a^) n \0) and 

I Qn+m 



(m\D(a)\n) = e 2 



VM.dz^dz. 



- <0|e 

n > 1 



(zi—a*)a Jz2+a)a i 



10) 



21=22=0 
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With (o|e (2l - Q * )a e (z2+Q)at |0) = e (* 1 - a *K»+ a ), in the case m > n the differentiation 
yields 



1 f) n 

(m\D(a)\n) = e^H 2 — _ — (z 2 + a)^"^ 

yrv.rm oz 2 



1.4.15) 



2 2 +a 

Comparison with the generating function of the Laguerre polynomials [1] yields 

n\ m _„ _ii„,i2 
— 7« 

ml 

(m\D^(a)\n)* = (-l) m - n (m\D(a)\n)*. (1.4.17) 




m\D(a)\n) = J — a m - n e—^ L™~ n (|a| 2 ) (1.4.16) 



A corresponding expression can be derived for n > m. 

1.4.3 Franck-Condon Factors in Electron Scattering and Tunneling 

We consider a single electron in one dimension, 

H = — + \(x)u (a + a ] ) + uj G a ] a. (1.4.18) 
2m 

Step change, \(x) = 9(x)\, A > 

In this case, the wave function for scattering states at energy E > can be written 

as 

oo 

*(x<0) = J2 ( a ^ lknX + bne~ iknX ) \n) 

n=0 

oo 

^(x>0) = ^t n e ik '- x \n) x (1.4.19) 

n=0 

with unshifted oscillator states |n)o in the expansion on the right side and shifted 
oscillator states \n)\ in the expansion on the left side of the step. Continuity of 
\^(x) and ^'(x) at x = yields two equations 

oo 

a m + b m = )J n (m\n)x 

n=0 

oo 

k m (a m -b m ) = ^k' n t n {m\n) x . (1.4.20) 



n=0 



We immediately recognize the appearance of the Franck-Condon Factors (m\n)\ 
that describe the sudden change of the oscillator potential when the electron moves 
across the step. 
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Adiabatic change of X(x) 

In this opposite case, X(x) changes very slowly as a function of x: let E > be 
the energy of a total scattering state (electron coming from the left). Then, the 
corresponding electron momentum is 



V 



V2mE, (1.4.21) 



and we write the time of flight tf through a typical distance Ax = X'(x)/X(x) over 
which the coupling potential X(x) varies as tf = mAx/p. Adiabaticity then means 
that the time of flight tf is large compared with 1/ujq, 

muin X'ix) 

> 1 



y/2mE X(x) 



-o » (1.4.22) 

Replacing E by the minimum total energy which is the oscillator ground state 
energy \TvjJq leads to 

2 Q Ax 2 (1.4.23) 

condition for the adiabatic situation. 
The adiabatic case is no longer simply describe by Franck-Condon factors. In 
contrast, the step change discussed above means Ax = across the step and is 
thus the opposite of the adiabatic situation Eq. (1.4.22). 



1.5 Franck-Condon Blockades and Noise for the 
U = oc-Anderson-Holstein Model 

1.5.1 Rate Equation Approximation 

The Rate Equation Approximation consists in neglecting all non-diagonal elements 
of the reduced system density operator. We thus replace 

(0,n\p(t)\0,n') -> 6 nn ,{0,n\p{t)\0,n) (1.5.1) 

in all terms in the Master equation. The validity or otherwise of such an approx- 
imation is in general hard to check - unless one directly compares results from 
a Master Equation and a Rate Equation calculation for a given model. Usually, 
one argues that for certain types of calculations off-diagonal elements might be 
small and may thus be neglected. Using the Rate Equation Approximation in the 
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Master Equation of the small Polaron model, we find a set of differential equations 
for the occupation probabilities 

Po, q (t) = (0,q\p(t)\0,q), P 1>q (t) = (l,q\p(t)\l,q), (1.5.2) 

where the label q is used here instead of n for the sake of a unified notation when 
comparing with recent results in the literature. The Rate Equations read 

p o, g = - E [ p o, ? rJ4, - P w r;c°J (1-5.3) 

g' 

K = - E fo^tS - ^r°4] , (1.5.4) 



where the transitions rates r°^L are defined as 



r££ = 2 \(q\X\q')\ 2 la [e + cu (q'-q)} (1.5.5) 

a=L,R 

r J4 = 2 E \{q\XW)\ 2 J a [uo(q' - q)) (1-5.6) 
7„(e) = m/ a (e)f a (e), %(e) = irv a (e)[l - f a (e)]. (1.5.7) 

The terms in Eq. (1.5.3) can be easily interpreted as follows: consider the level in 
the 'empty' state 0,q (no electron, phonon state q). This state changes by 



c di a i cxlc i 

in vibrational state q' . The energy difference 



- tunneling of an electron into the molecule at a rate rjf^L, leaving the molecule 



-Efinal - ^initial = £\q' ~ E 0q = E + U (q' - q) (1.5.8) 

has to be provided by the leads: before tunneling onto the molecule, the electron 
must therefore have had the single particle energy E& na \ — initial below the Fermi 
surface, which is why the rate is proportional to the occupation f a of the leads 
from which the electron tunnels. Energy conservation requires 

/affinal ~ ^initial) = fa + ^o(<?' ~ <?)) • (1.5.9) 

These processes lead to a decrease of the probability for the level being empty. 

- tunneling of an electron from the molecule in state 1, q' into the lead at a rate 
r^g, leaving the molecule in vibrational state q and thus increasing P q . The 
remaining energy, as given by the energy difference 

^initial - £fmal = £lq> ~ S 0q = S + U (q' - q), (1.5.10) 
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is now provided to the leads in the form of an additional electron with single 
particle energy initial — -^finai above the Fermi surface, which is why this rate is 
proportional to the Pauli-block factor 

1 " /^initial - ^final) = 1 - f a {g + LU (q' - q)) . (1.5.11) 

Note that in both cases the argument of the Fermi function is the same. Further- 
more, the tunneling density of states u a gives the number of possible single-particle 
states in the lead a with these energies. Finally, the squares of the Frank-Condon 
factors |(g|A|g')| 2 enter into all the transition rates: they describe the overlap of 
the vibrational wave functions according to the Frank-Condon principle. 

1.5.2 Current and Noise 

A typical current-voltage characteristic of the small polaron model with e = is 
shown in Fig. 1.5.2 a), which we took from the paper by Koch and v. Oppen [2]. 
For intermediate coupling strength A = 1, the current increases in steps ftwo when 
the transport window 

H L = V/2 > fJiR — -V/2 (1.5.12) 

is opened by increasing the source-drain (bias) voltage V across the molecule. For 
larger coupling strength A = 4, the current is strongly suppressed at small V. 
This is due to the exponential suppression of the Frank-Condon factors |(g|X|g')| 2 
at large A (shown as two-dimensional 'histograms' in Fig. 1.5.2 b) which means 
that only processes with a significant energy difference uo(q' ~ <?) contribute to 
transport, but at small V these are simply not allowed by energy conservation 
which therefore leads to a blockade of the current, the Frank- Condon blockade. 

Koch and v. Oppen found that Monte-Carlo (quantum trajectory) calcula- 
tions of the time-dependent rate equation Eq. (1.5.3) reveal deep insight into the 
transport process of the Anderson-Holstein model at strong coupling A = 4. At 
intermediate voltage bias V, the transport is strongly suppressed by the Frank- 
Condon factors which leads to long periods where no electrons pass through the 
molecule which remain in, say, state q = 0. At large A, the Frank-Condon factors 
exponentially grow with q for transitions with — > q. These rare transitions then 
trigger a chain of further transitions to still larger phonon numbers q, leading to 
an avalanche-type passage of a large number of electrons through the molecule in 
a very short space of time. These avalances are characterised by a hierarchical 
structure, i.e. when magnifying them as a function of time they display a similar 
structure on that new time-scale, again with long waiting times interrupted by 
'second generation' avalanches, and so on. This is shown in Fig. 1.5.2. 
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Figure 1.2: a) I-V curve for infinite-C/ Anderson-Holstein model at A = 1 and A = 4; 
b) Frank-Condon factors K^AT^')! 2 for A = 1 (left) and A = 4 (right); c) differential 
dl/dV plot at A = 4 and d) at A = 1. From Koch, v. Oppen, Phys. Rev. Lett. 94 
206804 (2005). 



1.6 Nonlinear Electron-Phonon Coupling 

We now discuss the case that is quite opposite to the linear electron-phonon cou- 
pling discussed so far. The material in this section is discussed in detail in a 
Diploma Thesis (Hannes Huebener, University of Hamburg 2007, unpublished) . 
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7 



time 



Figure 1.3: Avalanche type of transport for large A = 4: The number of electrons 
N(t) having passed through the molecule, and the phonon number q(t) have a self-similar 
structure as a function of time t. Parameters are T = 0.05u;o and V = 3ujq. From Koch, 
v. Oppen; Phys. Rev. Lett. 94, 206804 (2005). 

1.6.1 Model Hamiltonian 

Our model system is a resonator or a molecule that can be charged and un-charged 
by tunneling of single electrons. We assume that the resonator oscillates in re- 
direction with harmonic oscillator frequency ujq and mean square displacement Iq 
at zero temperature in the uncharged state. Charging with one additional electron 
leads to a Coulomb attraction towards a nearby metallic surface due to the image 
charge at distance 2a from the (uncharged) resonator equilibrium position, cf. 
Fig. 1.4. We assume strong electronic Coulomb blockade and thus two electronic 
states (|0), 'empty' and |1), occupied), leading to the effective Hamiltonian of the 
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Figure 1.4: (Left) Oscillator near a conductive surface with image charge effect. (Right) 
single electrons (de)charging the dot at rates 1^(1^) lead to a switching between a 
harmonic resonator potential (oscillator length Zq) and a non-linear potential due to 
the Coulomb attraction of the charged resonator towards the surface. From Huebener, 
Brandes; preprint (cond-mat). 



dot-resonator 

where M = H/IqUq is the effective oscillator mass of the resonator, u the unit 
vector from the dot to the metallic surface, and d) the creation operator of an 
electron on the dot, where for simplicity we disregard the electron's spin degree of 
freedom. 



1.6.2 Vibrational Eigenstates 

These have to be calculated by an exact diagonalization of H& ot . In the three- 
dimensional problem, the resonator motion in x-direction decouples from the mo- 
tion perpendicular to the x-axis for both charged and uncharged states, thus 
leaving us with an effectively one-dimensional problem: the uncharged vibra- 
tional states |0,n) simply are harmonic oscillator eigenstates with respect to the 
equilibrium position a. On the other hand, the charged states |l,n) are ob- 
tained by expansion into the eigenbasis of I = hydrogen atom eigenfunctions 

Mx) = ^ m h-i*xLU (f.) and subsequent numerical diagonalization. 

Here, x is measured from the metallic surface in units of a. 
The system Hamiltonian thus can be written as 

H s = ^ y]e i|W |w)(in|, \i, n) = \i) e ® \n) h (1.6.2) 

i=0,l n 

where \n)o are the harmonic oscillator eigenstates and |n)i are the vibrational 
eigenstates for the nonlinear potential V(x) in which the charged oscillator moves, 
cf. Fig. 1.4. 
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The model has two dimensionless parameters, 

n = = — e = e ' Ma = n A 

ft Zq ' 8ft 2 7re fto; ' 



where f2 determines the distance from the wall relative to the harmonic oscillator 
(zero-point fluctuation) length scale Iq, and e is given by Q times the typical 
Coulomb (charge - image charge) energy, 

E a = — e — , (1.6.4) 
Anemia 

relative to the harmonic oscillator phonon energy ftc^o. 

An important feature of the corresponding spectrum of the energy eigenvalues 
e vn for v = 1 is the appearance of an avoided level crossing (Fig. 1.5(a) inset 
near Q = 13.4). This can be traced back to the effective double-well structure of 
V(x) for coupling strengths e > 9. In this regime, the charged oscillator can be 
approximated by a two-level system at low energies, and one thus expects quantum 
coherent behavior due to tunneling between two states approximately localized 
around x = a and x = (i.e. near the surface). This vibrational tunneling can be 
made visible by electronic transport through the resonator, which requires the use 
of the full Master equations (and not the rate equations). 



1.6.3 Coherences and Master Equations 

As before, electrons jump off and onto the system according to 

H v = (V k J k J + H.c.), s=|0)(l|. (1.6.5) 

ka=L/R 

There is no dependence on the vibrational degree of freedom at this stage! This 
comes in at the point where we formulate the Master equation. As before, we use 
our general expression Eq. (1.3.17) and the procedure outlines thereafter: take 
matrix elements of Eq. (1.3.17) by 'sandwiching' with the exact eigenstates of Tis, 
insert the time-dependence of the s(t') etc., perform the Markov approximation in 
the interaction picture and transform back into the Schrodinger picture. The only 
difference is that now we do not use a polaron transformation but rather work with 
the exact eigenstates of the system Hamiltonian Tis- The result of this derivation 
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Figure 1.5: (color online), (a) Potential of an oscillator with image charge effect 
depending on the oscillator frequency and the eigenvalue spectrum of the correspon- 
dong Hamiltonian (for e = 12) (inset). Level anti-crossing occurs at £1 = 13.4. (b) 
Franck-Condon overlaps of the eigenfunctions of the charged and uncharged oscillator as 
a function of Q. There are two different regimes before and after the anti-crossing point. 

(c) dl/dV as a function of gate and bias voltage for e = 12 and f2 = 11 without the 
effect of the Franck-Condon factors {jij = 1) showing the usual conductance diamonds. 

(d) As (c) but with Franck-Condon factors: Current is surpressed for low bias and can- 
not be restored by changing the gate (Franck-Condon blockade). Negative differential 
conductance occurs as violet lines in the transport spectrum (c.f color scale) 



yields 



dt 



— H £ 0n' — £0n. 
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and a corresponding equation for the (l,n\p(t)\l,n'). Again, the Franck-Condon 
Factors 



are determined by the overlap of the vibrational wave functions before and after 
tunneling of one electron, but they now have to be calculated numerically. Some 
Franck-Condon factors are shown in Fig. 1.5(b). As a function of Q, there are two 
different regimes where either Tu or r 2 i is dominating. The sudden rise of T\\ is due 
to the level anti-crossing after which the shape of the first eigenstate resembles the 
harmonic oscillator ground state and vice versa for T21 and the second eigenstate. 
Also the sign of the has to be taken into account contrary to the usual rate 
equation approach where the Franck-Condon factors occur as squared values only. 

1.6.4 Vibrational Coherences and Electronic Noise 

The vibrational coherence (due to the level splitting) can be extracted from the 
electronic noise spectrum S(u). This quantity can also be calculated from the 
Master equation. We calculate the frequency dependent noise from the expression 
(which we do not derived in these Lectures) 



where Xr is the part of the Liouvillian that generates a quantum jump (an electron 
tunneling to the, e.g., right reservoir), TZ(u) is the pseudo-inverse of the Liouvillian 
£ and (•) = Tr[»p stat }. For simplicity, we only calculate the noise contribution in 
the right reservoir which (together with the contribution from the left reservoir) 
is sufficient to obtain the noise for strongly asymmetric capacitance coefficients. 

Figs. 1.6(b)&(c) demonstrate that the noise is indeed extremely sensitive to 
changes in the energy level spectrum of the resonator. The level splitting close 
to the avoided level crossing becomes visible as the appearance of a feature in 
S(u) at frequencies ui close to the corresponding Bohr frequencies (level splitting 
in the vibrational spectrum, Fig. 1.6(b) inset). The electronic noise thus acts as a 
direct probe of vibrational coherence. We point out that the appearance of non- 
trivial spectral information in S(u) is due to retaining the non-diagonal elements 
in the Master equation. In order to emphasize this point, we have compared our 
results with a corresponding calculation where these elements were set to zero (Fig. 
1.6(b)). (not discussed here). 

1.7 Electronic Damping of an Oscillator 

We now change our point of view and direct our attention on the oscillator. In 
our model with a single molecular level and a single phonon mode, we split the 



(0, n\l,m) 
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(a) (b) 




Figure 1.6: (color online), (a) Zero Frequency Fano factor for zero damping (lower 
frame) and several damping strength 7 G [10 -3 ,810 -3 ] (upper frame), (b) Frequency 
dependent Fano factor F(u) = S(u>)/(T) for systems with eigenfrequencies f2 as indicated 
in the inset, Vsias = Vbate = 30. Coupling parameters are Tl = 1/2tt and Tr = l/4ir. 
For comparison the case where only diagonal elements of the Liovillian are considered 
(black line), (c) As (b). The sensitivity of the noise to the level anti-crossing is seen 
as a growing and flattening peak at the frequency corresponding to the level spacing for 
eigenfrequencies scanning through the anti-crossing point. The oscillator eigenfrequency 
of each line is indicated by the right frame as well as its position with respect to the 
avoided-crossing. The scale of the axis refers to the oscillator frequency at the level anti- 
crossing point Q = 13.4 (red line) the other lines have an offset of ±n0.05. (d) F(uj) 
for £1 = 13.4 and several damping strength 7. For high damping the peak at the level 
spacing disappears and the noise spectrum resembles the single resonant level result. 



Hamiltonian into 'system' (oscillator) and 'bath' (electrons) according to 

H = Hs + Hsb + H b (1.7.1) 

Hs = u^a = ^ Fi + l -Mulq 2 (1.7.2) 

H B = H FA = ^2 £ k, a c{ a c ka + ^2(V ka c{ a d + H.c.) +ed ] d (1.7.3) 

k,a=L/R ka 

H S B = \'uj ( a + a ] )d ] d= \qd)d. (1-7.4) 
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Here, we neglect the electron spin in order to simplify things. We also renamed the 
coupling constant in 'a + at notation' in Tiss into A' so that we can use A as the 
coupling constant in q (oscillator position) notation. The 'bath' Hamiltonian 7ig 
is the Hamiltonian H.fa of the Fano- Anderson model, which is an exactly solvable 
model for electrons tunneling between two reservoirs a = L/R through a single 
level e. 

Rather than being interested in electronic transport in the presence of a phonon 
mode, we now focus on the oscillations of the phonon mode in presence of electron 
tunneling on and off the molecular level e. The idea is to 'integrate out' the 
electrons altogether and thus arrive at an effective theory for the oscillator. This 
programme can be carried out within the framework of Feynman- Vernon influence 
functional theory [3] 

1.7.1 Influence-functional Theory 

Let us consider the time-evolution of a density operator x{t) for a generic system- 
bath Hamiltonian, 

H = H S + H B + H SB , X(t) = e- iHt X (t = 0)e mt . (1.7.5) 

For the moment, let us assume that the system has one degree of freedom q and the 
bath the degree of freedom x (the generalisation to many bath degrees of freedom 
Xi is straightforward). We then use a representation of x(t) in spatial coordinates, 

{x,q\x{t)\q',x') = J dq dq' dx dx' (x,q\e~ lHt \q ,x )(x q \x(t = 0)\q' Q , x' ) 

x (x' q' \e im \q',x'). (1.7.6) 
We trace out the bath degree of freedoms to obtain an effective density matrix 

p(t) = Tr BX (t) (1.7.7) 

of the system, 

(q\p(t)W) = J dq dq' dx dx' dx{x,q\e~ lHt \q ,x )(x q \x( y t = 0)\q' , x' Q ) 

x (x' q' \e lHt \q',x). (1.7.8) 

Now we realise that the Hamiltonian H = Hs+Hb + Hsb induces a classical action 
S'totai = Ss[q] + Sb[x] + Ssb[xq], where in the following for notational simplicity 
we omit indices at the three S. We use the path integral representation for the 
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propagator matrix elements, 

(q\p(t)\q') = I dqodq'MpiOM) f exp [i (%] - S[q'])] F[q(t'), (/(?)] 

J Jq Jq' 

F[q(t'), q'(t')] = J dx dx' dx(x \p B \x f ) 

px PX 

x Vx V*x' exp [i (S[x] + S[xq\) - i {S[x'\ + S[x'q'])] 



•I'D 



• In the original Feynman- Vernon method, one assumes a factorising initial 
condition x(t = 0) = p(0) ® Pb, although that can be generalised to non- 
factorising initial density matrices x{t = 0); cf- f° r example H. Grabert, P. 
Schramm, G. L. Ingold, Phys. Rep. 168, 115 (1988), or the book by Weiss. 

• The functional q'(t')] is called influence functional . It describes 
the effect of the bath on the time-evolution of the system density matrix. 

• For zero system-bath coupling Hsb = 0, q'{t')} = 1 
Let us assume that we can write 

H SB = H SB [q} = f(q)X (1.7.9) 

with some given bath operator X and some given function f(q) of the system 
coordinate q. The influence functional can then be written as 

F[q(t'),q'(t')} = J dx dx' dx(x \p B \x' ) 

PX PX 

x / Vxexp[i(S[x} + S[xq])] / V*x' exp [-i (S[x'} + S[x'q'})] 

Jxo Jx' Q 



where U B [q] is the unitary time-evolution operator for the time-dependent Hamil- 
tonian H B + Hs B [q] with a given q(t'),0 < t' < t. Note that q(t') and q'{t') are 
independent paths, they enter as 'external' parameters into the influence func- 
tional which then in the final expression for (q\p(t)\q') is integrated over all paths 
q(t') and q'{t'). This form is useful to recognise general properties of jF[g(t'), q'(t')\, 

• q{t') = q'(t>)^F[q{t'),q'{t')] = l. 

. \T[q(t%q'(t')]\<l. 
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The 

Operator Form of Influence Functional 

FW)A'{t>)\ = Tr B (p B Ul[q']U B [q]) , (L711) 

is particularly useful for discussing the coupling to other baths (spin-baths, Fermi 
baths etc.) 



1.7.2 Non-equilibrium Damping: Quasi-Classical Approximation 

In the quasi-classical approximation, one expands around deviations from the 'di- 
agonal' in the double path integrals, i.e. paths q(t') = q'(t'). This is achieved by 
introducing center-of-mass and relative coordinates, 

x t = ^ (qt + q' t ) , yt = qt-q' t , (1.7.12) 

where we introduce a shorthand notation xt = x(t) etc. 
The system-bath coupling in our model then becomes 



HsB[qt>] = Ksb xv + -?r = \x t > + ^-Sd. (1.7.13) 



2 



Vt> 



2 



However, if the bath is the Fano- Anderson model, Tis = T~Cfa, something impor- 
tant is missing: the effect of the bath (i.e. the electrons) on the oscillator will 
depend on the oscillator coordinate q, because the electronic energy depends on 
q as (e + Xq)Sd. This would mean that we actually have to regard the coupling 
term 7i.sB[qt'} as part of the bath. The 'bath' 

HpA + Xqd^d, (1.7.14) 

however, is a very difficult bath as the effective level energy wildy fluctuates, and 
there is no hope to integrate out the electrons with Eq. (1.7.14) exactly. 

Mozyrsky, Hastings and Martin [4] therefore came up with the idea to use an 
adiabatic approximation (Born-Oppenheimer approximation) for the case where 
the oscillator moves much slower than the electrons tunnel on and off the level, 

u Q <^T R ,r L ,V, (1.7.15) 

where V is the bias voltage. 

Let us follow the time-evolution of the oscillator over a short time interval, 
starting at its density matrix center-of-mass coordinate Xq at time t = 0. Let us 
assume we just want a very short time interval [0, t] in which the oscillator barely 
moves, in which however many electronic tunnel events happen. The electrons will 
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permanently try to shift the oscillator equilibrium position via the polaronic shift, 
in particular if the coupling A is strong. We then write 



Hb + 7~Lsb [Qf 



Ufa + Xqt'd^d = Ufa + Xxod'd + A (t'x + 



W 
2 



dU 



Hi 



Hsi 



(1.7.16) 

which means that we shift the 'main part' xq of the oscillator coordinate into the 
bath Hamiltonian 7Yb which thus acquires the desired dependence on the oscillator 
coordinate, and we shift the 'fluctuating' part (deviations from xq and the whole 
quantum off-diagonal coordinate y t > into the system-bath coupling Hsb- 

We will now do a 'self-consistent' calculation: first, calculate the double path 
integral and the influence functional with fixed xq in 7Yb; then in the resulting 
expressions replace the fixed Xo again by the full c.o.m. coordinate Xt- This 
is a form of Born-Oppenheimer approximation: we first keep the 'nuclei' (the 
oscillator) fixed ('clamped nuclei approximation') and are thus able to calculate 
the nuclear potential generated by the electrons, and then use the parametric 
nuclear coordinates as full nuclear coordinates in the so generated potential. 

The propagator J for the density matrix p of the oscillator defines the time- 
evolution, 



J(x,y,t;x ,y ) 



dx dy (x + — \p(0)\x 



yo 

2 



J(x,y,t;x ,q ) 



Vx / Vyexp 



dt' (Mx t >y t > - Mulxfyt') 

x F[x v ,y v ] (1.7.17) 
The influence functional in a second order cluster expansion in A is 

InFlxt^y?] = -iX / dt'y t '(n(t')} B 

Jo 



fdt' f ds {/+ - f t 7} {S(f - 8 )f+ - S*(t' - a)fr] 
Jo Jo 



ft 



A 



t ; x 



... ± 



S(t) = (n d (T)n d (0)) 



'\.7AX 



The only additional approximation now is the Born-Oppenheimer approximation, 
i.e. we neglect the terms .... in f t , = A [t'cto + ... ± . Then the influence 
functional becomes 

\nF[x t >,yt>] = 



X 2 [ dt' [ dsy t ,y s ReS(t' - s) 
Jo Jo 

2iX 2 x I dt' I dsy t ,sImS(t' - s). (1.7.19) 



-iX j dt'y t >(n d ) B 
-t' 
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We will now reckognize the physical meaning of these expressions: we can re- 
arrange the whole propagator into 



J BC (x,y,t; x ,2/o) 



■to 



Vx Vyexp [iM(x t y - x y )] x 



x exp 



\ dt'y t ,K(t') -X 2 dt! / dsy t ,y s Re(h d (t' - s)h d (0)) 
Jo Jo Jo 



B 



K{t') = Mx t > + Mu 2 Q x t > + X(h d ) b + 2 [ ds\ 2 slm(h d (t' - s)n d (0)} § x 

Jo 

(1.7.20) 

This is an interesting expression: it is essentially a 'huge' Gaussian integral cen- 
tered around K = with a width described by quadratic, i.e. Gaussian fluc- 
tuations, the strength of which is determined by the real part of the correlation 
function (n(t ; — s)n(0))^. The condition K — now is the desired differential 
equation of our oscillator, 

Mx t > + Muj 2 x t , + X(h d ) s + 2 f ds\ 2 slm(n d (t' - s)n d (0)) s x = 0, (1.7.21) 

Jo 

where — A(n) ^ is an extra force on top of Mu^Xf. We can calculate the occupation, 
e.g. with the help of Green's functions, 



(n d )(x) 



dw T L f L (u) + T R f R (uj) 



2tt (u-e- Ax) 2 + \(T L + T R y 
For T L = T R = T, zero temperature T = and fiL = V/2 = —fi R , we have 

du 6(V/2 + 6(-V/2 - u) 



1.7.22) 



1 

"2k 



2tt (u-e- \x) 2 + T 2 

e + Xx - V/2 e + Xx + V/2 
arctan h arctan n 



1.7.23) 

1.7.24) 
1.7.25) 



In the spirit of the Born-Oppenheimer approximation, one would now replace the 
constant xq and xo by time-dependent expressions, thus leading to 

Mx t > + Mu 2 x t > + X(n d ){x t ') + 2 [ dsX 2 slm(n d (t' - s)n d (0)} B x s = 0. (1.7.26) 

Jo 

We see that there is then a retarded friction force coupling to the velocity x s of 
the oscillator. The oscillator furthermore moves in an effective potential, U e g(x) 
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and is subject to an effective force F eff = —U' eS {x) 

U cS (x) = \Mulx 2 + [ dx'X(n d }(x'). (1.7.27) 
* Jo 

Depending on the position of the level e, the bias voltage V, the coupling strength 
A, the tunnel rate T, and the bare mass and frequency of the oscillator, this 
potential can give rise to interesting bistabilities. 



2. MANY PHONON MODES AND 
DISSIPATION 



Bulk phonons are often a source of dissipation when they interact with electrons. In 
metals, the electron-phonon coupling is usually weak. Electron-phonon scattering 
rates approach zero when the temperature T goes to zero. This is because at T = 0, 
no occupied electronic states states are available above the Fermi level Ep and all 
states below E F are occupied. Electron-phonon scattering would be only possible 
by absorption of phonons which however is proportional to the Bose distribution 
ub{u)) = [exp(hw/kBT) — which vanishes (u > 0) at zero temperature. 

In electronic few-level systems, electron-phonon scattering is however possible 
even at zero temperature in the form of spontaneous emission of phonons. In 
electron transport experiments, this has been observed, for example in coupled 
quantum dots. 

2.1 Phonons in Double Quantum Dots 

Double quantum dots are simple charge qubits and thus candidates for quantum 
logic gates. Let us construct a coupled quantum dot model by taking into account: 
1. strong interactions 2. coupling between the dots 3. coupling to the external 
world (leads) 4. coupling to other degrees of freedom such as bosonic excitations, 
photons, phonons. 

No-one has succeeded to reliably calculate transport properties of such a model 
in the whole parameter space. In the following, we again use our Master equation 
formalism and treat the coupling to the leads in lowest order perturbation theory 
(T ^> Tjc, no Kondo physics). Coulomb effects are also dealt with in an extremely 
simplified manner: Assume two dots coupled by a tunnel barrier. The Coulomb 
charging energy U in small dots is typically U ~ 1 meV 3> transport = (jU£ — A t i?)/ e: 
assuming there are N electrons in the left and M electrons in the right dot 
ground state |0) = \N,M). The additional transport electron is either on the left 
or on the right dot ~* states \L) — \N + l,M),\R) — \N,M + 1). Therefore, only 
three states are involved! Nevertheless, including dissipation, microwave radiation 
etc this still is complicated enough to be an interesting model. 
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2.1.1 Isolated Double Dot 

We describe the space spanned by the two states L and R with the help of (pseudo) 
spin Pauli matrices: The Hamiltonian for this subspace is 

H = ^a z + T c a x , (2.1.1) 

where a z = \L)(L\ - \R){R\ and a x = \L)(R\ + \R)(L\. Here, e = e L - e R is the 
difference of the ground state energies of \L) and \R), and T c describes the coupling 
between the two dots. 

The eigenvectors and energy eigenvalues of H are 

|±> = -^-[±2T c \L) + (A T e)\R)}, N ± := a/ 4|T c | 2 + (A =f e) 2 

l\l± 

e± := ±^A, A:= y/ e* + 4\T C \* . (2.1.2) 
In microwave experiments, the square-root form of A was detected in 1998. 

2.1.2 Coupling to Phonons 

We describe the coupling to phonons by a Hamiltonian 

Q Q 

Here, ujq are the frequencies of phonons, and the gQ denote interaction constants. 
Although not exactly solvable, the model is quite well understood for closed sys- 
tems [5] (isolated dots with one additional electron). The coupling to external 
leads offers the possibility to study its non-equilibrium properties, such as the 
inelastic stationary current through the dots. 

2.1.3 Equations of Motion 

We describe the dynamics of the double dot by a reduced (effective) statistical 
operator p(t) . The coupling to the reservoirs is described by the additional Hamil- 
tonian 

H rcs + H v , H ies = J2 £ k c l c k + T, £ % d l d k ( 2 - L4 ) 

k k 

H v = ( V k c k s L + W k d\s R + c.c.) , s L :=\0){L\, s R :=\0)(R\. 

k 

(2.1.5) 
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We start from the Liouville-von Neumann equation for the density operator of the 
total system (DOT + PHONONS + RESERVOIRS), 

± x (t) = -i[H, X (t)}. (2.1.6) 

The first step: transformation into the interaction picture with respect to H : = 
Hb + H rcs , 0(t) = e tH ^ t Oe~ %H ° t for any operator O. Write the total Hamiltonian 
H = H + V, and Eq. (2.1.6) becomes 

j t m = -mt),xm (2-1-7) 

Iteration of Eq. (2.1.7), 



j t m = -mt),x(o)]- 



[2.1.8) 



We define the effective density operator for the double dot, 

pit) = Tr reSjPhonon x(t). (2.1.9) 

Now, second order Born and Markov approximation in the perturbation V yields 
[6,7] 

d 

-q^PLL it) = -iT c \p LR (t) - PR L (t)]+T L [l- p L L(t) ~ pRR{t)\ 

d 

g^PRR(t) = ~lT c \p RL {t) - PLR{t)\ ~ FRPRR{t), 

where the tunnel rates are 

T L = 27r^|y fc | 2 5(4-£ i ) r R = 27rJ2\Wk\ 2 S(e L k -e R ). (2.1.10) 

k k 

and we allow for tunneling from a left reservoir at rate Tl into the left dot, and 
from the right dot to the right reservoir at rate T R . 

For the remaining equation for the off-diagonal element p^ R = p* RL , one has 
to choose between a perturbation theory in g Q (weak coupling, PER), or a per- 
turbation theory in T c in a polaron-transformed frame (strong coupling, POL) . 
In general, no exact solution of the model is available: this is the case even for 
coupling to one bosonic mode only (g<g oc 5q,q , Rabi Hamiltonian). 

For the spin-boson problem with T r /l = 0, it is well-known that POL is 
equivalent to a double-path integral 'non-interacting blip approximation' (NIBA) 
that works well for zero bias e = but for e ^ does not coincide with PER at 
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small couplings and low temperatures. It turns out that when comparing both 
approaches for T R / L ^ 0, one finds nearly perfect agreement for very large e T c , 
a regime that has been tested experimentally recently [8] . 

The standard Born and Markov approximation with respect to A yields 

j/l^) = [ie-^ P -r R /2}p LR (t) + [iT c -S_}p RR (t)-[iT c -5 + }p LL (t). 
Here, the rates are 

7p := 2vr^p( A) coth (/3A/2) , 5 ± := -§»p(A) coth (/3A/2) =F §£p(A) 

p(u,) := ^> Q | 2 %-u, Q ). (2.1.11) 
Q 

where A := e 2 + 4T C 2 is the energy difference of the hybridized levels, and (3 = 
l/k R T the inverse phonon equilibrium bath temperature. Note that beside the 
off-diagonal decoherence rate j p , there appear terms oc 5± in the diagonals which 
below turn out to be important for the stationary current. 

On the other hand, the polaron transformation [6] leads to an integral equation 



2 c(t - t') PLR (t') 



p P L ° R L (t) = - I dt'e^ 
Jo 

+ iT c (t') {C(t - t')p LL (t') -C*(t- t')p RR (t')} 

where 



C(r) :=exp|-^ C 



duo^- [(1 - cos cut) coth(/5c/2) + % sinu;t] > . (2.1.12) 
us 1 



2.1.4 Stationary Current 

Calculated as / sta t = —e2T c lmpL R (z = 0) from Laplace transforming the equations 
of motion; 

t - -eT L T R G+ (o W\\ 

stat " T L G- + (T L + T R )G + -T L T R 

has to be used with either the perturbative result (small e-p coupling) or the 
polaron transformation (small inter-dot coupling), 

Gi PER) := 2r c Im , lTc ~ % - (2.1.14) 
te - 7 P - T R /2 
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or 



(POL) 



2T c Im 



-iT c C e 



G (POL) . = 2TcIm 



-iT c C% 



#.1.15) 



+ 



i + (i/2)r fl c e ' 



1 + (l/2)T R C y e 



where C £ := J °° dte l£t C(t). 
2.1.5 Noise 

The Master equation formalism has also been used to calculate the current noise 
spectrum S(u) and the Fano Factor F = S(0)/2eI. Recent experiments on 
transport through layers of double quantum dots seem to be in good agreement 
with these calculations (preprint G. KieBlich et al.). 
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